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Abstract 

We test if the latest Gold set of 182 SNIa or the combined "Platinum" set of 192 SNIa from 
the ESSENCE and Gold sets, in conjunction with the CMB shift parameter show a preference 
between the ACDM model, three luCDM models, and the DGP model of modified gravity as 
an explanation for the current accelerating phase of the universe's expansion. We consider flat 
wCDM models with an equation of state w{a) that is (i) constant with scale factor a, (ii) varies 
as Wo + Wa(l — a) for redshifts probed by supernovae but is fixed at —1 at earlier epochs and 
(iii) varies as uiq + Wa(l — a) since recombination. We find that all five models explain the data 
with comparable success. 



Observations of Type la supernovae (SNIa), the cosmic microwave background, and large scale 
structure corroborate that the expansion of our universe is accelerating. Explanations of this 
acceleration often invoke the existence of a dark energy component with the unusual property of 
negative pressure [1]. The simplest example of dark energy is the cosmological constant A with an 
equation of state w = p/ p = —1. Although in the context of A, p and p are the effective pressure 
and energy density attributed to vacuum energy, in general they are the corresponding quantities of 
the fluid describing dark energy. Possibilities with a time-varying equation of state are modelled by 
introducing a scalar field called quintessence \2\. No dark energy model has a natural explanation. 
Explanations of the current acceleration that avoid the introduction of dark energy emerge from 
modifications to gravity, typically in the context of braneworld models. 

For models with dark energy, we consider the nonflat ACDM model, and three flat universe 
models with equations of state that differ from — 1 at some time in the expansion history. We refer 
to these as flat wCDM models. The simplest of these has a constant w that differs from —1. To 
model dark energy with a time- varying w, we adopt a parameterization that is well-behaved at 
high redshift: w{a) = wq + Wa{^ — a) [2]) where a is normalized such that it is equal to 1 today and 
is related to redshift zhj a = 1/(1 -|- z). While this parameterization reduces to the linear relation 
'w{z) = WQ+WaZ [4| at low redshifts, it has the disadvantage that when analyzing high-redshift data, 
we are implicitly restricting ourselves to models that are well-represented by this time evolution. 
Said differently, we are placing the strict prior that dark energy must have evolved in this manner 
throughout. Keeping the latter in mind, one case we consider is that dark energy evolves according 
to w{a) = wq + Wa{^ — a) for redshifts probed by SNIa and behaves like a cosmological constant at 
earlier epochs. A second case, involving a stronger prior, is that w{a) = wo + Wa{^ — a) is valid till 
the epoch of last scattering. 

The Friedman equation in dimensionless form with H = a/a and Hq = H{z = 0) is 

= n„(l + z)3+0^(l + ^)3(l+-o+-.)e-3»..T^+j^^(l + ^)2^ 

where ilrm and 0^ are the matter, dark energy and curvature densities in units of the critical den- 
sity. For the ACDM model, fi^ = wq = —1 and Wa = 0. For flat models, 0^ = 1 — — = 0. 

The DGP model [5] (named so for its authors) is a generally-covariant infrared modification of 
general relativity (GR) and is not reducible to an extension of GR with additional scalar or vector 
degrees of freedom. At short distances gravity is 4-dimensional, while in the far infrared, gravity 
appears to be 5-dimensional because of gravitational leakage from the 4-dimensional brane into the 
5-dimensional bulk. The implications for cosmology are that at early times the correction from the 
infrared modification is negligible and the universe obeys the standard cosmology. However, at late 
times, the weakening of gravity is significant and leads to self-accelerated expansion without the 
need for any form of matter [6] . Due its firm theoretical foundation |[7] , we analyze the DGP model 
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as the canonical example of modified GR that explains the current acceleration without the need 
for dark energy. 

The Hubble expansion is given by [B] 

= [yWr + + + Z)3]2 + 17,(1 + zf . (2) 

Here Q,r = i^^c^o)^^ ^ where Vc = M|/(2M|) is the length scale beyond which 4-dimensional 
gravity (with Planck scale M4) transits to 5-dimensional gravity (with Planck scale M5). Setting 
z = in Eq. yields fi^ = 1 — [^/Tl^. + ^/Qr + ^m]^- Note that the DGP model has the same 
number of parameters as ACDM, with fi^ replacing 

To compare the various models on an equal-footing, we only analyze data that probe the ex- 
pansion history of the universe. We do not utilize data that are sensitive to the evolution of density 
perturbations since these have not been fully worked out for the DGP model, although progress has 
been made in Ref. [8]. We analyze the distance moduli of the Gold set of 182 SN [9] (of which 16 
have z > 1) compiled from Refs. \10\ \TT\ [12] and the "Platinum" set of 192 SN [TS^ compiled from 
the ESSENCE [14j and Gold sets. In our analyses, we include the CMB shift parameter [15] which 
measures the distance to the last scattering surface. For most of what follows, we do not use the 
baryon acoustic oscillation (BAO) distance parameter A extracted from the scale corresponding to 
the first acoustic peak at recombination [T6]. The procedure used to determine A assumes that w 
does not vary with redshift. The resulting value may not be applicable for time- varying w [17] . 
making it unsuitable for three of the five models we are considering. Recent joint analyses of the 
DGP model with older data in various combinations (with some including the A parameter) can 
be found in Ref. [IE]. 

The shift parameter, defined in terms of the i?o-independent luminosity distance = HqcIl 
(where di \s the luminosity distance) [15] : 

R^V^ f,f''''\ , (3) 

(1 + ZCMB) 

is approximately equivalent to the ratio of the sound horizon at recombination to the comoving 
distance to the last scattering surface. For all practical purposes, R is model- independent. We 
use the value R = 1.70 ± 0.03 [E] obtained from the WMAP 3-year data [20] with a redshift at 
recombination zcmb = 1089 it 1. 

The statistical significance of a model is determined by evaluating Xr = {R"^'^ — R*^)"^ /cr"^ and 
xItv' which after marginalization over a nuisance parameter, has the absolute value [21] 

xIn = A-^. (4) 
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where 

Here, and ai are the distance modulus and its uncertainty at redshift Zi. N = 182 for the Gold 
set and = 192 for the Platinum set. 

ACDM: In Fig.[T]we display the results of our analysis of the nonflat ACDM model. The shaded 
regions are the Icr, 2(t and Scr allowed regions from an analysis of the latest SN data (upper panel: 
Gold set, lower panel: Platinum set) and of the CMB shift parameter R. Although expected, the 
orthogonality of the regions is striking. The solid contours depict the corresponding regions from 
the joint analysis. The best-fit parameters are provided in Table 1. The dot marks the best-fit from 
the joint analysis. Notice that the combined analyses prefer universes that are almost flat with a 
tendency for positive curvature. If we restrict ourselves to flat universes, the minimum value 
for the joint analysis with the Gold set (Platinum set) increases to 163 (196.3), and the best-fit 
moves to $7^ = 0.3 {Q.m = 0.25); xIat increases to 158.6 for the Gold set but remains essentially 
unchanged for the Platinum set. 

Flat lyCDM with constant w: We allow w to take values different from —1, but do not allow 
for time variation by setting = 0. Figure [2] is similar to Fig. [1] except that we plot wq vs Vtm- 
Keeping in mind that the number of degrees of freedom is the same as that for the ACDM model, 
we see that the fit to the ACDM model is not significantly better; with the Gold set (Platinum set) 
it has a minimum that is only 1.8 (0.3) lower. 

Flat twCDM with late-varying w: In this case, w has time-variation only from z = 1.8 until 
today, and equals —1 at earlier times. We also require 0.15 < Vtm < 0.35. In Fig. [21 we only show 
the la and 2a regions for the separate SN and R analyses for obvious reasons. It is no surprise 
that the R parameter shows no sensitivity to Wa because we have assumed that the dark energy 
behaves as a cosmological constant since recombination until z = 1.8. This assumption eliminates 
the constraining power of R that comes from its long lever arm. The joint constraint is dominated 
by the SN data. This is evident from Fig. [3] and Table 1. The allowed regions and best-fit point do 
not move significantly on adding the R parameter to the analysis. Although the Gold and Platinum 
sets favor different regions of Wa, on comparing the allowed regions with the regions occupied by 
different dark energy models in the {wq^Wo) plane, as classified in Ref. [22j, it is clear that no class 
of models is excluded even at the 2a C. L. by either dataset. 

Flat wCDM with varying w: The equation of state is allowed to vary since recombination. 
Again, we require 0.15 < ilm ^ 0.35. As expected, the results of analyzing SN data alone are 
identical with that of late-varying w since there are no SN data with z > 1.8. See Table 1 and 
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Gold 


(182 SN) 




Platinum 


(192 SN) 




ACDM 




Qa 




x' 


Qa 




SN 


156.4 


0.95 


0.48 


195.2 


0.85 


0.33 


SN+i? 


158.4 


0.68 


0.36 


195.6 


0.74 


0.27 


SN+i?+BAO 


161.4 


0.72 


0.30 


195.6 


0.74 


0.27 


u;CDM, Wa = 


x' 


Wo 




x' 


Wo 




SN 


156.6 


-1.75 


0.46 


195.4 


-1.16 


0.31 


SN+i? 


160.2 


-0.85 


0.28 


195.9 


-0.94 


0.24 


SN+i?+BAO 


160.3 


-0.86 


0.29 


196.6 


-0.98 


0.26 


wCBM, w{z > 1.8) = -1 


x' 


Wo 


Wa 


x' 


Wo 


Wa 


SN 


156.5 


-1.11 


2.39 


195.3 


-1.11 


-1.16 


SN+i? 


156.5 


-1.28 


2.69 


195.5 


-1.06 


0.81 


u;CDM 


x' 


Wo 


Wa 


x' 


Wo 


Wa 


SN 


156.5 


-1.11 


2.39 


195.3 


-1.11 


-1.16 


SN+i? 


157.0 


-1.35 


1.54 


195.5 


-1.09 


0.67 


DGP 


x' 






x' 






SN 


156.4 


0.24 


0.36 


195.1 


0.22 


0.24 


SN+i? 


160.3 


0.14 


0.23 


196.4 


0.16 


0.17 



Table 1: The minimum values and best-fit parameters for each of the five models. The psr 
degree of freedom is substantially different for the Gold and Platinum datasets. All tt;CDM models 
are flat. In the analyses in which both wq and Wa are allowed to vary freely, we rec[uire that 
take values between 0.15 and 0.35. All five models have Xr = at the minimum in the analysis 
of the shift parameter alone. We do not show the corresponding best-fit parameters because the 
x'ji distributions are too broad for the parameters to be meaningful. We have included the BAO 
constraint for the two models with a constant w for comparison. The number of parameters in the 
analyses of the successive models are 3, 3, 4, 4, 3, respectively, where in each case one parameter is 
a nuisance parameter that is marginalized over. 

Figs. [3] and m From Fig. 01 we now see how R helps to constrain Wa- Nevertheless, all the dark 
energy models classified in Ref. [22] remain safe. With respect to ACDM, the additional free 
parameter improves the minimum x'^ of the joint analysis with the Gold set (Platinum set) by only 
1.4 (0.1). 

DGP: Figure [5] is similar to Fig. [1] with JIa replaced by the physical parameter relevant to 
DGP, ^Ir- A similar orthogonal relationship exists between the SN and R regions. However, here 
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we see that the joint analysis with the Gold set (Platinum set) prefers a slightly open universe with 
= 0.027 ±0.014 {Qk = 0.041 ±0.010). The overall fit is only slightly worse compared to ACDM. 



To assess the impact of the BAO constraint on the two models that do not have a time-varying 
equation of state (ACDM and the u)CDM model with constant w), we use the measured value 
A = 0.469 ± 0.017 dS], where 



to calculate Xa = {^"^^ — A*^)^/cj^. Here, zbao = 0.35 is the typical redshift of the SDSS sample 
of luminuous red galaxies. From Figs. [6] and [71 it can be seen that the BAO constraint provides 
a satisfying confirmation that the different datasets are concordant, and helps to further constrain 
the regions from the joint analysis of SN data and R. The minimum and best-fit parameters 
from the joint analyses including the BAO constraint are provided in Table 1. For the ACDM 
model, the addition of this datapoint to the joint analysis with the Gold set results in a preference 
for a universe with less curvature at the expense of increasing the minimum by 3. The BAO 
constraint has no effect on the best-fit parameters of the ACDM model in a joint analysis with the 
Platinum set. 

In conclusion, it is noteworthy that all 5 models fit the SN data equally well. Only on inclu- 
sion of R in the analyses do minor differences develop. Current data cannot tell if the accelerated 
expansion is caused by a cosmological constant, by dark energy with a constant w, by dark energy 
whose equation of state started varying recently or has always been varying, or due to modified 
gravitational physics of conventional mattei0. A detailed understanding of how density perturba- 
tions evolve in braneworld cosmologies will enable the use of the vast amount of data available on 
the power spectrum from observations of the CMB and large scale structure and may help with the 
basic question of whether the acceleration is due to a new form of energy or due to new aspects of 
gravity. A step forward from kinematical probes of modified gravity to dynamical ones is the order 
of the day. 
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^ Other analyses of similar datasets can be found in Ref. [23] . 
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ACDM 




Figure 1: The la, 2a and 3a allowed regions for the ACDM model from an analysis of SN data 
(upper panel: Gold set, lower panel: Platinum set) of the CMB shift parameter R and from a joint 
analysis. The best-fit parameters from the joint analyses are indicated by a dot. See Table 1. 
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Figure 2: Similar to Fig. [H but for the wCDM model with a constant equation of state w{a) = wq. 
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Flat w'CDM with m^(z>1.8)= -1 

8\ 




Figure 3: The allowed regions for the flat u;CDM model with an equation of state w{a) = wq + 
Wa{l — a) in the redshift range < 2 < 1.8 and a constant equation of state w = —1 ior z > 1.8. 
Only the Icr and 2a regions are shown for separate SN and R analyses, while the 3a region is also 
shown for the joint analysis. The best-fit parameters from the joint analysis are indicated by a dot. 



Flat vt^CDM 




Figure 4: Similar to Fig. [Sj but with an equation of state w{a) = wq + Wail — a) for < z < 1089. 
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Figure 5: Similar to Fig. [H but for the DGP model. 
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Figure 6: Similar to Fig. [H but also showing the effect of the BAO constraint on the parameter 
space. 
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Figure 7: Similar to Fig. [21 but also showing the effect of the BAO constraint on the parameter 
space. 



15 



